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Abstract 



We calibrate and test various variants of field theory models of the interest rate with data from 
eurodollars futures. Models based on psychological factors are seen to provide the best fit to the 
market. We make a model independent determination of the volatility function of the forward rates 
from market data. 



1 Introduction 

In this paper, we compare field theory models of interest rate models with market data, and propose 
certain modified models inspired from theoretical considerations and observed facts about the interest 
rates. The theoretical framework for all these models is Baaquie's formulation [0, || of forward rates 
as a two dimensional quantum field theory. The Baaquie model is a generalization of the Heath- Jarrow- 
Morton (HJM) model; the key feature of the field theory model is that the forward rates f(t, x) are 
imperfectly correlated in the maturity direction x > t, and which is specified by a rigidity parameter //. 
The models we study are the following: (a) forward rates with constant rigidity (b) forward rates 
with the variation of the spot rate constrained by a new parameter S , and two new models proposed in 
this paper, namely (c) forward rates with maturity dependent rigidity [l(x — t), and lastly (d) forward 
rates with non-trivial dependence on maturity specified by an aribitrary function z = z(x — t). 

We first briefly review Baaquie's field theory model and review the market data used in this study. 
We then test the field theory model, introduce two variants and test them as well. We find that the 
correlation structure can be explained by a relatively straightforward two parameter model which also 
has a useful theoretical interpretation. 



2 



The HJM model 



2.1 



Definition of the model 



In the HJM- model the forward rates are given by 




(1) 



where Wi are independent Wiener processes. We can also write this as 




(2) 



where rji represent independent white noises. The action functional, is 

K 



s[w\ = -\Y j J dtnm (3) 



We can use this action to calculate the generating functional which is 



e 4 EE* ( 4 ) 



3 Field theory model with constant rigidity 

We now review Baaquie's field theory model presented in Q with constant rigidity. Baaquie proposed 
that the forward rates being driven by white noise processes in (^) be replaced by considering the forward 
rates itself to be a quantum field. To simplify notation, we write the evolution equation in terms of the 
velocity quantum field A(t,x), and which yields 

f(t,x) = f{t ,x) + / dt , a(t',x) + Y, dt'a l (t',x)A l (t',x) ih) 

J to i— i J to 



or 

K 



d.f(t,x) 
dt 



= a(t, x) + <Tj(t, x)Aj(t, x) (6) 



The main extension to HJM is that A depends on x as well as t unlike W which only depends on t. 

While we can put in many fields Ai, we will see that the extra generality brought into the process 
due to the extra argument x will make one field sufficient. Hence, in future, we will drop the subscript 
for A. 

Baaquie further proposed that the field A has the free (Gaussian) free field action functional 

with Neumann boundary conditions imposed at x = t and x = t+TpR. This makes the action equivalent 
(after an integration by parts where the surface term vanishes) to 

1 f°° ft + TFR / 1 Q 2 \ 

s= ~2j t dt ] t dxA ^ x )y---^Q^) A ^ x ) ( 8 ) 



This action has the partition function 



Z[j] = exp (j^ 



ti rt+T FR \ 

dt dxdx'j(t,x)D(x -t,x' -t)j(t,x')\ (9) 



with 

cosh /i (TpR -\B- 6'\) + cosh/x (T FR -(6 + 6')) 
1 ' ' FR) ^ 2smh^iT FR (10) 

= D(9', 6; T FR ) : Symmetric Function of 9, 9' 

where 9 — x — t and 9' = x' — t. We can calculate expectations and correlations using this partition 
function. Note that due to the boundary conditions imposed, the inverse of the differential operator 
D actually depends only the difference x — t. The above action represents a Gaussian random field 
with covariance structure D. In Q], a different form was found as the boundary conditions used were 



Dirichlet with the endpoints integrated over. This boundary condition is in fact equivalent to the 
Neumann condition which leads to the much simpler propagator above. In the limit Tfr — ► oo which 
we will usually take, the propagator takes the simple form /ie~ M0> cosh ^9 < where 9 > and 9 < stand for 
max(0, 9') and mm(8, 9') respectively. 

When fi — > 0, this model should go over to the HJM model. This is indeed seen to be the case as 
it is seen that lim^o D(6, 9'; Tfr) — 5^-. The extra factor of Trr is irrelevant as it is due to the 
freedom we have in scaling <r and D. The er we use for the different models are only comparable after 
D is normalized^]. On normalization, the propagator for both the HJM model and field theory model in 
the limit ji — ► is one showing that the two models are equivalent in this limit. 

The basic model with constant rigidity can be generalized in many different ways. The generalization 
to positive valued forward rates, and to models with stochastic volatility are studied in ||. In this paper 
we generalize the free field model to more complex dependence of /j, and f(t, x) on the maturity direction 
9. 



4 The Market Data used for the Study 

We used the Eurodollar futures data for the following study. A Eurodollars futures contract is represents 
a deposit of US$1,000,000 for three months at some time in the future. Currently, futures contracts for 
deposits upto ten years into the future are actively traded. Significant historical data for contracts on 
deposits upto seven years into the future are available. If one makes the reasonable approximation that 
f(t, 9) is linear for 9 between contract times, one can use this data as a direct measure of the forward 
rates. 

Further, the straightforward simplification that the Eurodollar futures prices directly reflect the for- 
ward rate was done, an assumption previously used in the literature [Q. This is a reasonable assumption 
as the forward rates are small enough that the difference between the logarithmic measure of the forward 
rate used in theory and the arithmetic rates used in the market are insignificant. We also attempted to 
analyse Treasury bond tick data from the GovPx database but we found it impossible to obtain forward 
rates accurate enough for our purposes. The main reason for this is that while we were able to obtain 
reasonably accurate yields for a few maturities, the differentiation required to get the forward rates 
from the yields introduced too many inaccuracies. This is somewhat unfortunate since Treasury bonds 
represent risk free instruments while a small credit risk exists for Eurodollar deposits. 

For the following analysis, we used the closing prices for the Eurodollar futures contracts in the 1990s. 
This is exactly the same data as used by Bouchaud Q as well as Baaquie and Srikant || . In Bouchaud 
0, the spread of the forward rates and the eigenfuctions of its changes in time are analyzed. For our 
purposes, we found it more useful to look at the scaled multivariate cumulants of the changes in forward 
rates for different maturity times. 



5 Assumptions behind the tests of the models 

The main assumption that has to be made for all the tests of the models is that of time translation 
invariance. In other words, we have to assume that a(t, 9) is actually only dependent on 9 and not 
explicitly on t. We also assume that the propagator D(9,9') has no explicit time dependence which is 
possible in principle. It is reasonable and conceptually economical to assume that different times in the 
future are equivalent. Further, carrying out any meaningful analysis while these quantities are subject 
to changes in time is impossible. 

Another important assumption that has to be made is that the forward rate curve is reasonably 
smooth at small intervals at any given point in time. This assumption is very difficult to test in any 
meaningful sense given the relative paucity of data as forward rate data is available only at 3 month 
intervals (which is what necessitates this assumption in the first place). However, the assumption is a 
reasonable one to make as one would intuitively expect that the forward rate, say three years into the 

1 This freedom exists since we can always make the transformation u(8) ~ ri(9)a(0) and D(8,9') ~ D(Q,0')/(r](6)r](8')) 
without affecting any result 




Figure 1: The lines of constant 9 for which we have obtained the forward rates by linear interpolation 
from the actual forward rates which are specified at constant x. 



future would not be too different from that three years and one month into the future. In fact, we will 
show later that there seems to be strong evidence of very long term correlations in the movements of the 
forward rate. This seems to make the smoothness assumption reasonable as nearby forward rates tend 
to move together (except possibly at points very close to the current time). This assumption is required 
as the forward rate data is provided for constant maturity which we have been denoting by x while we 
want data for constant 9, as shown in figure 0. With this assumption, we can get the data by simple 
linear interpolation. The loss in accuracy due to this linear interpolation is not all that serious if e, the 
time interval of t between specifications of the forward rates is small as the random changes which we are 
interested in will be much larger than the introduced errors. This same procedure was used in Matacz 
and Bouchaud j| as well as Baaquie and Srikant ||. 
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Figure 2: The empirically determined function a{9). 



6 The Correlation Structure of the Forward Rates 

A very interesting quantity to look at in the analysis of forward rates f(t, 9) is the correlation (or scaled 
covariance) among their changes for different 9. Specifically we are interested in the correlation between 
Sf(t, 9) and 8 fit, 9'), where Sf(t, 9) = f(t + e, 9) — f(t, 9). Using a free (Gaussian) quantum field theory 
model, this quantity should be equal to 

c(e Q/) = {Sf(t,9)Sf(t,9'))-{Sf(t,9)){Sf(t,9')) = D{9,9') 

y/{6p(t,0)) - (5f{t,9))^{5p{t,9')) - y/D(9,9)D(9\9') 



To a reasonable degree of accuracy, we can ignore the first order expectations such as (8f(t, 9)) as they are 
much smaller than the second order expectations if e is small. For an e of one day, the error is completely 



negligible especially given the other approximations. We will do so for the rest of the chapter. If we have 
a model for the propagator D{9,9'), we have a prediction for this correlation structure. Alternatively, 
we can use the correlation structure to fit free parameters in D{9, 9'). 

It should be noted that for free (Gaussian) quantum fields the correlation is independent of cr(#), so 
no assumption of its form has to be made. This is the reason why we used the scaled covariance rather 
than the covariance itself to perform the study. It is equivalent to fixing the inherent freedom in the 
quantities a and D to make D(8,8) = 1. The reduction in the freedom of a also allows us to directly 
estimate it from data since we have <j(8) = \J < 5f 2 (t, 8) > if D(8,8) = 1. This is shown in figure |[ 
Further, the correlation between innovations in the forward curve is given exactly by D. The correlation 
structure in the market estimated from the Eurodollar futures data is shown in figure 0. The structure 
is fairly stable in the sense that the correlation structure for different sections of the data are reasonably 
similar. 

Since the propagator is always symmetric, it will be convenient to calculate only D(9 < ,9 > ) for the 
different models where 9 < = mm(9,9') and 9 > — max(0, #') for purposes of comparison. 
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Figure 3: The correlation structure observed in the market. 
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Figure 4: Observed correlation structure : a different view. 

For the one factor HJM model, this correlation structure is constant as all the changes in the forward 
rates are perfectly correlated. In other words, D(9 7 8') = 1. For the two factor HJM model, the predicted 



Covariance of the changes in the forward rates 




Figure 5: The covariance of innovations of forward rates observed in the market 



correlation structure is given by 

c(g 0I) = vimvim+vtfw) = i+g(%(g ; ) (12) 

We see that this correlation structure depends on a function of g(6) = °* Yl . Hence, a whole function 
has to be fitted from the correlation structure, something which is quite infeasible. The covariance might 
be a better quantity to test the two factor HJM model as the prediction of the covariance has a simpler 
form 



C(9,9') = a 1 (e)ai(d') + a 2 (9)(T2(e') (13) 

We still need to specify a functional form for <j\ and 02 as it is not possible to estimate entire functions 
from data. The usual specification of <Ji(9) = <7o and <J2{9) = a\e~ xe inspired by the assumption that the 
spot rate follows a Markov process is easily seen to be unable to explain many features of the covariance 
in figure |^ such as the peak at one year or the sharp reduction in the covariance as the maturity goes 
to zero. We can straightaway conclude that the one factor HJM model is insufficient to characterize the 
data while the two factor HJM model provides us with too much freedom as we can put in an entire 
arbitrary function to explain the correlation structure. If we try to reduce the freedom by theoretical 
considerations, we are again unable to explain the data. 

We will see that the field theory model with constant rigidity, while explaining some features of the 
correlation, does not predict the correlation very well. Hence, we consider generalizations to the constant 
rigidity model. 



7 Analysis of Field Theory Model with Constant Rigidity 

We have analysed this model in detail in the previous section. We have seen that the model describes 
the innovations in the forward rates in terms of a Gaussian random field A whose structure is defined by 
the action in (^). For convenience, we repeat the action below in terms of the variables t and 9 = x — t 

(14) 




To obtain the predicted correlation structure from the propagator (|Io|), we have to take the limit 
TpR — > oo and obtain 

D{9, 9') = cosh f i6 < = | (e-^ 9 -^ + e -^ e+e '^ (15) 

The predicted correlation structure for this model can be found from this form of the propagator by 
normalization and from (|ll|) is given by 



^W^w (16) 

when the limit TpR — > oo is taken. To estimate the parameter \i from market data, we use the Levenberg- 
Marquardt method from Press et al [|| to fit the parameters to the observed correlation structure graphed 
in figure ||. The fitting was done by minimizing the square of the error. The overall correlation was fitted 
by fi = 0.061/year. To obtain the error bounds, the data was split into 346 data sets of 500 contiguous 
days of data each and the estimation done for each of the sets. The 90% confidence interval for this 
data set is (0.057, 0.075). Note that the confidence interval is asymmetric from the overall best fit due to 
the nonlinear dependence of the correlation (|l^) on /i. The root mean square for the correlation for the 
best fit value is 4.23% which shows that the model's prediction for the correlation structure is not very 
good. The main problem as can be seen from a comparison between the prediction for the best fit [i in 
figure ^ and the actuall correlation structure in figure [| is that the prediction is largely independent of 
the actual value of 9 and largely determined by \9 — B'\ which is not the case in reality. The correlation 
rapidly increases as 9 increases in reality. 

Filted Correlation for Constant Rigidity Model 
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Figure 6: Fitted correlation for constant rigidity model 



Filled Correlation lor Constrained Constant Rigidity Model 
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Figure 7: Fitted correlation for constrained field theory model 



8 Field Theory with Constrained Spot Rate 

One clear fact we notice from the covariance of the innovations in the forward rates in figure || is that the 
covariance falls rapidly as 9 — ► 0. This observation leads one to a model g where A(t, 0) is constrained 



to follow a normal distribution with variance a. The mean of A(t,0) can be fixed at any value but will 
cause a corresponding change in a(0) which makes the mean value irrelevant. For calculational purposes 
it is easiest to assume that it remains at zero. This constraint can be implemented by modification of 
the action to 



d^Se 1 ' 



where S is the action specified in ([14|). The propagator D(9, 9') for this model is given by 



D(9,9') = \ie 



cosh fi8 < — 



fie 



H + a 



After normalizing, we see that the prediction for the correlation structure is given by 

0(0,9') 



■ (cosh fi8 < 



^ e -M0<(cosh M 0>-H£^) 



(17) 



(18) 



(19) 



We can see that the free parameters are n and a. Further, it will be seen that it is easier to consider 
the ratio a//j, 2 as it is dimcnsionless. The results of the Levenberg-Marquardt method showed that the 
fitted value of \i and a were very small, of the order of 10 _7 /year for /i and 10 _13 /year 2 for a both 
being very unstable but the ratio a//i 2 was stable with a value in the range (6.7, 10.7) with an overall 
best fit of 9.4. The most reasonable explanation for this behaviour is that the ratio a/// 2 determines the 
behaviour of ( |l9| ) for small /i and it is this region of parameter space that gives a correlation structure 
closest to the empirically observed one. We see from the fitted propagator in figure that the behaviour 
at large 9 is slightly better when the constraint is put in. The root mean square error was 3.35% which 
again means the fit was not very good though significantly better than if the constraint was not applied. 
It must be recognized that the constraint introduces one extra free parameter which should improve 
the best fit. Hence, we see that this model, while again performing better than HJM, is still not very 
accurate. While the results are not very good, they do represent a reasonable first approximation and 
are still significantly better than the one factor HJM model. 



9 Field Theory Model with Maturity Dependent Rigidity fi = 

Another way to get a correlation structure that depends directly on the values of 9 and 9' in a significant 
way and not only on their difference is to make fi a function of 9. This has a direct physical meaning 
as it means that if we imagine the forward rate curve as a string, its rigidity is increasing as maturity 
increases making the A for larger 9 more strongly correlated if /i decreases as a function of 9. We 
choose the function /j, = jt^§ as it declines to zero as 9 becomes large as is expected from the observed 
covariance in figure |5[ contains the constant n case as a limit and is solvable. The action is given by 

This is still a quadratic action and can be put into a quadratic form by performing integration by parts 
and setting the boundary term to zero since we are assuming Neumann boundary conditions. The inverse 
(Greens function) of the quadratic operator or the propagator for this action is found to be 

D(9,9';T FR )- M ° a 



2Xa(a + 1/2)(1 - (1 + XT FR )- 2a ) 

+ AT ™)~ 2a ( 1 + Xd >) a ~ 1/2 + (i + A#>r Q - i/2 ) (21) 
(Sy^ 1 + XTF ^ 2a{ - 1 + ^<r~ 1/2 + (i + am— 1 / 2 ) 



where a = y j + and where we have put the bound on the 9 variable Trr explicitly. The reason for 
this is that the limits have to be taken carefully in order to compare this model to the HJM in the limit 
/io — > and to the constant rigidity field theory model when A — > 0. 

Let us first consider the limit A — > 0. First, we note 



1 Mo Mo L , A 2 mo 
4 + A 2 ~ Ay 4/i 2 ~ A 



«=t/7 + ^-T^ + 73~T (22) 



Therefore, we have 

(1 + Afl)-"- 1 / 2 = ((1 + A(9) 1/A ) ^° (1 + A<9)~ 1/2 ~ e-^° 9 (23) 

Similarly (1 + A<9) Q ~ 1/2 - e^° e , (1 + A6>)~ Q ~ 1/2 - e^ Q0 and (l + \T Fn )- 2a ~ e - 2w,T ™. Putting all these 
limits into ([H]) and performing some straightforward simplifications, we see that (|2l|) becomes equal ( |ic| ) 
in the limit A — > 0. In the taking of this limit, we did not have any trouble with Tfr. However, for the 
HJM limit, we will see that the limit Tfr — * oo has to be taken only after the limit /io — > has been 
taken. 

Let us now consider the limit /io - * 0. In this limit a ~ ^ + jS-- Hence, only one term in ( pl| ) survive 
as all the others are multiplied by a — 1/2. This surviving term can be evaluated 

2A a -1/2 l-(l + AT Ffl )-i (1 + ATi?fl) 

= m| x 2A^ x 1 + XTpR 1 

2A /i 2 ATf_r 1 + XTfr 
1 



Tfr. 



The terms (1 + A#>) Q-1 / 2 and (1 + A6*<) Q_1 / 2 obviously go to one in this limit and so were not included in 
the calculation above. This result can be seen to be equivalent to the HJM propagator after normalization. 
If the limit Tfr — > oo is taken first, then the propagator becomes 

which exhibits a 9 dependence in the limit /io — > 0. Hence, this cannot be made equivalent to HJM if 
the limits are taken in the wrong order. This problem is not present in the constant rigidity model. 

For comparison with market data, we still take the limit Trr — > 00 as the model is then still directly 
related to the field theory model. The predicted correlation structure for this model is then given by 

C(9 9>\ A" + y2)(l + Afl < ) 2 °+q-l/2 \ * 



We fitted the parameters /zo and A to the correlation structure observed in the market in a similar 
manner as for the field theory model and obtained the results /io = 1.2 X 10 _5 /year and A = 0.108/year. 
The root mean square error in the correlation was 3.35%. On performing the error analysis for the 
parameters, it is found that /io is very unstable but always very small (less than 10 _2 /year) while the 
90% confidence interval for A is (0.099, 0.149). The relatively high value for A seems to show that the 
falloff of the rigidity paramater fi — j-^g is fairly rapid. The error is reduced from 4.23% to 3.35% but 
an extra parameter has had to be added and the model has become considerably more complicated due 
to the freedom of the form of the rigidity parameter /i. Further, we seem to be in the region of very 
small fiQ which does not behave well in the HJM limit. In fact, the correlation structure in this limit is 
given by 



(27) 



Fitted Correlation for Non-constant Rigidity Model 




Figure 8: Fitted correlation structure for the non-constant rigidity model. 

Due to the very small value of /i for the fitted function, this is a very good approximation for the fit. 
The obtained fit for the correlation function can be seen in figure |^. 

The limited improvement, the relatively complicated form of the correlation and the near zero /xo 
problem prompted us to consider a different way of approaching the problem which presented a much 
more satisfactory solution. This model is described in the next section 

10 Field Theory Model with f(t, z(6)) 

To see where we might make an improvement, we notice that the predicted correlation structure with the 
field theory model is largely defined by the e~^ e> ~ e< ^ term which means that the correlation does not 
depend explicitly on the times 9 > and #<0. However, we see immediately from figure [| that the correlation 
increases significantly as we increase #> and 9 > . This is intuitively reasonable as market participants 
are likely to treat the difference between ten and fifteen years into the future quite differently from the 
difference between now and five years. Far out into the future, we would expect all times to be equivalent. 
In other words, there is good reason to expect lime < _ >00 D(6 > , 0<) = 10. This is not satisfied by the 
constant rigidity models or by the varying rigidity model (if the limit Tfr — > oo is taken). For the latter 
model, this is slightly surprising since y, — > as 9 — > oo and we might expect that for large 9 the varying 
rigidity model should go into the HJM model limit (D = 1). However, this does not happen as previously 
discussed since we have taken the limit Tfr — > oo. 

Further, the relatively marginal reduction of the error shows that varying the rigidity parameter does 
not quite reflect the data. An alternative way to consider the problem would be to use the observed 
correlation structure to induce a metric onto the 9 direction. In some sense, this metric would be 
measuring the "psychological distance" in the investor's minds which corresponds to a certain separation 
in maturity time. To make this concrete, let us write the observed correlation as D{9,9') — e~ s ( e,e \ 
Since D(9, 9) = 1, s(9, 9) = and s is symmetric as well. If we can show the triange property (which in 
the case of one dimension reduces to the straightforward condition that s(9i, 63) = s(#i, #2) + s(6*2, #3)), 
we can see that s makes a good definition of distance in 9. From the market data, it can be shown that 
this rule is very approximately satisfied and we can use it as an approximate way to induce a metric 
onto the 9 direction from the observed market data. 

It should also be noted that introducing the metric is different from changing the form of the rigidity 

2 There is another term of the form e~ ' J ( 9 >+ e <) but this has only a small effect on the correlation structure 
3 Obviously, 9 < — > 00 automatically implies 0> — » 00 



function n(9). To sec this, we write the action with the rigidity function ji(6) as 



*--5jC*r*(^ + ?(s)j <»> 

where the functional variation of /x with has been absorbed into the variable z = g{9) (where g is 
invertible) so that the \x above is a constant. With a change of variables we get the action as 

s °"=-J,. dt L i;W r?W ) <29) 

where h — g . With the introduction of the metric, we obtain the action 

i fti /-ff(oo) / 1 / (94 \ 2 \ 

S neW = -o / dt dz U 2 + -2 «- ^ S °« (30) 



The Green's functions for S new should be solved using the z variables, and as expected the solution 
is given by D(z, z') — ^(exp — n\z — z'\ + cxp — fi(z + z 1 )). It can be shown that the martingale condition 
is satisfied with the Green's function given by D(z, z'). 

Bearing in mind the condition that, at large 9, the correlations should be close to 1, or equivalently 
that the distance should be small, we choose a metric that satisfies property, g(9) = tanh (36. We use 
this form of the metric to fit the correlation structure and obtain the result that /i = 0.48/year and 
f3 = 0.32/year with a root mean square error of only 2.46%. Both parameters are also stable when the 
error analysis for the parameters is carried out. The 90% confidence interval for [i is (0.45, 0.58) and 
that for f3 is (0.22, 0.33). Hence, we see that even the parameter estimation for this model is more robust 
as the parameters are atleast stable. Further, the shape of the fitted function is clearly closer to the 
observed one as can be seen from figures |[ [?], || and ||. The error that remains is largely confined to 
the correlation between the spot rate and other forward rates which is not too surprising since the spot 
rate behaves very differently from the other forward rates. 

Fitted Correlation tor metric model 
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Figure 9: Fitted correlation for the model with metric g(9) = tanh/30 

We emphasize here that this involves a fundamentally new way of thinking of the interest rate models. 
So far, we have made models which generalized HJM so as to achieve a theory without too little freedom 
as in the one factor HJM model or too much freedom as in the two factor HJM model. While retaining 
this framework, we now use empirical data to guide us in refining the model to give us an insight into 
market psychology which will result from the induced metric. 
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